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Abstract — Tradeoff in diversity, multiplexing, and delay in 
multihop MIMO relay networks with ARQ is studied, where the 
random delay is caused by queueing and ARQ retransmission. 
This leads to an optimal ARQ allocation problem with per- 
hop delay or end-to-end delay constraint. The optimal ARQ 
allocation has to trade off between the ARQ error that the 
receiver fails to decode in the allocated maximum ARQ rounds 
and the packet loss due to queueing delay. These two probability 
of errors are characterized using the diversity-multiplexing- 
delay tradeoff (DMDT) (without queueing) and the tail proba- 
bility of random delay derived using large deviation techniques, 
respectively. Then the optimal ARQ allocation problem can be 
formulated as a convex optimization problem. We show that the 
optimal ARQ allocation should balance each link performance 
as well avoid significant queue delay, which is also demonstrated 
by numerical examples. 

I. Introduction 

'in a multihop relaying system, each terminal receives 
the signal only from the previous terminal in the route 
and, hence, the relays are used for coverage extension. 
Multiple input-multiple output (MIMO) systems can provide 
increased data rates by creating multiple parallel channels and 
increasing diversity by robustness against channel variations. 
Another degree of freedom can be introduced by an automatic 
repeat request (ARQ) protocol for retransmissions. With the 
multihop ARQ protocol, the receiver at each hop feeds back 
to the transmitter a one-bit indicator on whether the message 
can be decoded or not. In case of a failure the transmitter 
sends additional parity bits until either successful reception 
or message expiration. The ARQ protocol provides improved 
reliability but also causes transmission delay of packets. Here 
we study a multihop MIMO relay system using the ARQ 
protocol. Our goal is to characterize the tradeoff in speed 
versus reliability for this system. 

The rate and reliability tradeoff for the point-to-point 
MIMO system, captured by the diversity-multiplexing trade- 
off (DMT), was introduced in [I]. Considering delay as the 
third dimension in this asymptotic analysis with infinite SNR, 
the diversity-multiplexing-delay tradeoff (DMDT) analysis 
for a point-to-point MIMO system with ARQ is studied in 
[2], and the DMDT curve is shown to be the scaled version 
of the corresponding DMT curve without ARQ. The DMDT 
in relay networks has received a lot of attention as well (see, 
e.g., [3].) In our recent work [4], we extended the point-to- 
point DMDT analysis to multihop MIMO systems with ARQ 
and proposed an ARQ protocol that achieves the optimal 
DMDT. 
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The DMDT analysis assumes asymptotically infinite SNR. 
However, in the more realistic scenario of finite SNR, re- 
transmission is not a negligible event and hence the queueing 
delay has to be brought into the picture (see discussions in 
[ ]). With finite SNR and queueing delay, the DMDT will be 
different from that under the infinite SNR assumption. The 
DMDT with queueing delay is studied in [5] and an optimal 
ARQ adapted to the instantaneous queue state for the point- 
to-point MIMO system is presented therein. 

In this work, we extend the study [5] of optimal ARQ 
assuming high but finite SNR and queueing delay in point- 
to-point MIMO systems to multihop MIMO networks. This 
work is also an extention our previous results in [4] to 
incorporate queueing delay. We use the same metric as that 
used in [5], which captures the probability of error caused by 
both ARQ error, and the packet loss due to queueing delay. 
The ARQ error is characterized by information outage prob- 
ability, which can be found through a diversity-multiplexing- 
delay tradeoff analysis [2], [4]. The packet loss is given 
by the limiting probability of the event that packet delay 
exceeds a deadline. Unlike the standard queuing models for 
networks (e.g., [6], [7]) where only the number of messages 
awaiting transmission is studied, here we also need to study 
the amount of time a message has to wait in the queue 
of each node. Our approach is slightly different from [5], 
where the optimal ARQ decision is adapted per packet; we 
study the queues after they enter the stable condition, and 
hence we use the stationary probability of a packet missing 
a deadline. An immediate tradeoff in the choice of ARQ 
round is: the larger the number of ARQ attempts we used 
for a link, the higher the diversity and multiplexing gain we 
can achieve, meaning a lower ARQ error. However, this is at 
a price of more packet missing deadline. Our goal is to find 
an optimal ARQ allocation that balances these two conflicting 
goals and equalizes performance of each hop to minimizes 
the probability of error. 

The remainder of this paper is organized as follows. 
Section II introduces system models and the ARQ protocol. 
Section III presents our formulation and main results. Nu- 
merical examples are shown in Section IV. Finally Section 
V concludes the paper. 

II. Models and Background 

A. Channel and ARQ Protocol Models 

Consider a multihop MIMO network consisting of N 
nodes: with the source corresponding to i = 1, the destination 
corresponding to i = N, and i = 2, ■ • • ,N—l corresponding 
to the intermediate relays, as shown in Fig. 1. Each node 



is equipped with Mj antennas. The packets enter the net- 
work from the source node, and exit from the destination 
node, forming an open queue. The network uses a multihop 
automatic repeat request (ARQ) protocol for retransmission. 
With the multihop ARQ protocol, in each hop, the receiver 
feeds back to the transmitter a one-bit indicator about whether 
the message can be decoded or not. In case of a failure 
the transmitter retransmits. Each channel block for the same 
message is called an ARQ round. We consider the fixed 
ARQ allocation, where each link i has a maximum of ARQ 
rounds i = 1, • • ■ N — 1. The packet is discarded once 
the maximum round has been reached. The total number of 
ARQ rounds is limited to L: ^^J^ 1 — ^- This fixed ARQ 
protocol has been studied in our recent paper [4]. 



(iv) 



We consider both the long-term static channel, where 
Hn = Hi for all I, i.e. the channel state remains 
constant during all the ARQ rounds, and independent for 
different i. Our results can be extended to the the short- 
term static channel using the DMDT analysis given in 
[4]. 
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Fig. 2: Left: full duplex multihop relay network. Right: half duplex 
relay multihop MIMO relay network. 
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Fig. 1: Upper: relay network with direct link from source to 
destination. Lower: multihop MIMO relay network without direct 
link. 

Assume the packets are delay sensitive: the end-to-end 
transmission delay cannot exceed k. One strategy to achieve 
this goal is to set a deadline ki for each link i with 
YliLi^ ki — Once a packet delays more than ki it 
is removed from the queue. This per-hop delay constraint 
corresponds to the finite buffer at each node. Another strategy 
is to allow large per-hop delay while imposing an end-to-end 
delay constraint. Other assumptions we have made for the 
channel models are 
(i) The channel between the ith and (i + l)th nodes is given 
by: 



SNR 



HjiXj 



W, 



l<l<L t . (1) 



The message is encoded by a space-time encoder into 
a sequence of L matrices G C MiXT ,: I = 

1, • ■ • , L}, where T is the block length, and Yij G 



C 



M i + 1 xT 



1, • • • , N — 1, is the received signal at 



the (i + l)th node, in the Zth ARQ round. The rate 
of the space-time code is R. Channels are assumed 
to be frequency non-selective, block Rayleigh fading 
and independent of each other, i.e., the entries of the 
channel matrices Hij G C Mi+lxMi are independent 
and identically distributed (i.i.d.) complex Gaussian with 
zero mean and unit variance. The additive noise terms 
Wi t l are also i.i.d. complex Gaussian with zero mean 
and unit variance. The forward links and ARQ feedback 
links only exist between neighboring nodes. 

(ii) We consider both the full-duplex and half-duplex relays 
(see, e.g., [4]) where the relays can or cannot transmit 
and receive at the same time, respectively, as shown 
in Fig. 2. Assume the relays use a decode-and-forward 
protocol (see, e.g., [4]). 

(iii) We assume a short-term power constraint at each node 
for each block code. Hence we do not consider power 
control. 




Fig. 3: The logarithm of the cost function (7) for the (4, 1, 2) 
multihop MIMO relay networks. SNR is 20 dB. 




Fig. 4: Allocations of optimal ARQ: L\, L%, LI+L2, kf, in a (4, 1, 
2) MIMO relay network. SNR is 20 dB. (The optimal fc| = k-kl.) 



B. Queueing Network Model 

We use an M/M/l queue tandem to model the multihop 
relay networks. The packets arrive at the source as a Poisson 
process with mean interarrival time fx, (i.e., the time between 
the arrival of the 71th packet and (n — l)th packet.) The 
random service time depends on the channel state and is 
upper bounded by the maximum ARQ rounds allocated Li. 
As an approximation we assume the random service time at 
Node i for each message is i.i.d. with exponential distribution 
and mean L;. With this assumption we can treat each node as 
an M/M/l queue. This approximation makes the problem 
tractable and characterizes the qualitative behavior of MIMO 



multihop relay network. Node i has a finite buffer size. The 
packets enter into the buffer and are first-come-first-served 
(FCFS). Assume p > Li so that the queues are stable, i.e., 
the waiting time at a node does not go to infinity as time 
goes on. Burke's theorem (see, e.g., [7]) says that the packets 
depart from the source and arrive at each relay as a Poisson 
process with rate Pi/pj, where pi is the probability that a 
packet can reach the ith node. With high SNR, the packet 
reaches the subsequent relays with high probability: pi ps 1 
(the probability of a packet dropping is small because it uses 
up the maximum ARQ round.) Hence all nodes have packets 
arrive as a Poisson process with mean inter-arrival time p. 

C. Throughput 

Denote by b the size of the information messages in bits, 
B[t] the number of bits removed from transmission buffer 
at the source at time slot t. Define a renewal event as the 
event that the transmitted message leaves the source and 
eventually is received by the destination node possibly after 
one or more ARQ retransmissions. We assume that under 
full-duplex relays the transmitter cannot send a new message 
until the previous message has been decoded by the relay at 
which point the relay can begin transmission over the next 
hop (Fig. 2a.) Under half duplex relays we assume transmitter 
cannot send a new message until the relay to the next hop 
completes its transmission (Fig. 2b.) 

The number of bits B transmitted in each renewal event, 
for full-duplexing B = (N — 1)6, and for half-duplexing 
B = (N— 1)6/2 when N is odd, and B = Nb/2 when N is 
even. The long-term average throughput of the ARQ protocol 
is defined as the transmitted bits per channel use (PCU) [2], 
which can be found using renewal theory [8]: 
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where r is the average duration from the time a packet arrives 
at the source until it reaches the destination node, and = 
denotes asymptotic equality. A similar argument as in [2] 
shows that E(t) = (N - 1)T for high SNR. 

D. Diversity-Multiplexing-Delay Tradeoff 



We cannot assume infinite SNR because otherwise the queue- 
ing delay will be zero, as pointed out in [5]. However we 
assume high SNR to use the DMDT results in our subsequent 
analysis. 

III. Diversity, multiplexing, and delay tradeoff 

VIA OPTIMAL ARQ ROUND ALLOCATION 
A. Full-Duplex Relay in Multihop Relay Network 

1) Per-Hop Delay Constraint: The probability of error 
depends on the ARQ window length allocation Li, deadline 
constraint ki, multiplexing rate r, and SNR p. For a given r 
and p, we have 

P e ({L i },{k i }\p,r) = 

N-l 

P ARQ { P , {Lj}) + ]T P Qucuc (A > h). (4) 

i=l 

Here Di denotes the random delay at the ith link when the 
queue is stationary. This P e expression is similar to that 
given by Equation (33) of [5]. Our goal is to allocate per-hop 
ARQ round {Li} and delay constraint {k^} to minimize the 
probability of error P e . 

For the long-term static channel, using the DMDT analysis 
results [4] we have: 
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Here fi(r) is the diversity-multiplexing tradeoff (DMT) for 
a point-to-point MIMO system formed by nodes i and i + 1. 
Assuming sufficient long block lengths, fi(r) is given by 
Theorem 2 in [1] quoted in the following: 

Theorem 1: [I] For sufficiently long block lengths, the 
diversity-multiplexing tradeoff (DMT) f(r) for a MIMO 
system with Alt transmit and M r receive antennas is given 
by the piece-wise linear function connecting the points 
(r, (M t - r){M r - r)), for r = 0, • • • , min(M t , M r ). 

Denote the amount of time spent in the ith node by 
the ?ith message as D l n . The probability of packet loss 
Poueuc{Di > ki) can be found as the limiting distribution of 
lim^oc P(D l n > ki) (adapted from Theorem 7.4.1 of [8]): 

Lemma 2: The limiting distribution of the event that the 
delay at node i exceeds its deadline ki, for M/M/l queue 
models, is given by: 



iWA > h) = lim P{D\ > k % 
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The probability of error P e in the transmission has two 
sources: from the ARQ error: the packet is dropped because 
the receiver fails to decode the message within the allocated 
number of ARQ rounds, denoted as P ARQ , and the probability 
that a message misses its deadline at any node due to large 
queueing delay, denoted as P Que u C . We will give P e for various 
ARQ relay networks. Following the framework of [1], we P e {{Li\,{ki\\p,r) = 
assume the size of information messages b(p) depends on the j = i 
operating signal-to-noise ratio (SNR) p, and a family of space 
time codes {C p } with block rate R(p) = b(p)/T = rlogp. 
We use the effective ARQ multiplexing gain and the ARQ 
diversity gain [2] 
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Here the difference in the service rate and packet arrival rate 
- > and utility factor — both indicate how "busy" 



node m is. Using the above results, (4) can be written as 
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Note that the queueing delay message loss error probability is 
decreasing in Lj, and the ARQ error probability is increasing 
in Li. Hence an optimal ARQ rounds allocation at each node 
Li should trade off these two terms. Also, the optimal ARQ 
allocation should also equalize the performance of each link, 
as the weakest link determines the system performance [ ]. 



Hence the optimal ARQ allocation can be formulated as 
the following optimization problem: 



„ min P e ({L 4 },{fc;}|p,r) 



where 



A = { 1 <Li <n, i = 1,. 
Si=i fet < fe- 
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The following lemma (proof omitted due to the space limit) 
shows that the total transmission distortion function (15) is 
convex in the interior of A. 

Lemma 3: The transmission distortion function (15) is 
convex jointly in Li and ki in the convex set 



{Li}, {ki} : ki > — jj. 



i = 1.---N - 1. 



Lemma 3 says that except for the "corners" of A the cost 
function is convex. However these "corners" have higher 
probability of error: ki and Li take extreme values and hence 
one link may have a longer queueing delay then the others. 
So we only need to search the interior of A where the cost 
function is convex. 

To gain some insights into where the optimal solution 
resides in the feasible domain for the above problem, we 
present a marginal cost interpretation. Note that the proba- 
bility of error can be decomposed as a sum of probability 
of error on the zth link. The optimal ARQ rounds allocated 
on this link should equalize the "marginal cost" of the ARQ 
error and the packet loss due to queueing delay. For node i, 
with fixed ki, the marginal costs (partial differentials) of the 
ARQ error probability, and the packet loss probability due to 
queueing delay, with respect to Li are given by 
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>0. (11) 



Note that f- < 0. The optimal solution equalizes these 
two marginal costs by choosing Li <E Note that 

these marginal cost functions are monotone in Li, hence the 
equalizing L* exists and 1 < L* < p if the following two 
conditions are true for Li = 1 and Li = p: 



(i): 



(ii) 
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These conditions involve nonlinear inequalities involving 
p, p, r, Mi and M 1+1 , which defines the case when the 
optimal solution is in the interior of A. Analyzing these 
conditions reveals that these conditions tend to satisfy at 
lower multiplexing gain r, small Mi or Afj+i, small ki, and 
larger p (light traffic). Note that with high SNR condition 



(ii) is always true for moderate k values. When (i) and (ii) 
are violated, which means one error dominates the other, 
then the optimal solution lies at the boundary of A. With 
the total ARQ rounds constraint in (8), using the Lagrangian 
multiplier an argument similar to above still holds. 

2) End-to-End Delay constraint: When the buffer per 
node is large enough a per hop delay constraint is not needed, 
and we can instead impose an end-to-end delay constraint. 
The exact expression for the tail probability of the end-to- 
end delay is intractable. However a large deviation result is 
available. The following theorem can be derived using the 
main theorem in [9]: 

Theorem 4: For a stationary M/M/l queue tandem (with 
full-duplex relays): 



lim lim — log P ( 
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where 9*= mm?^ 1 



This theorem says that the bottleneck of the queueing net- 
work is the link with longest mean service time L t . Hence 
the optimal ARQ round allocation problem can be formulated 

as: 



min P e ({Li},{ki}\p,r) 



(14) 
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For high SNR, this can be shown to be a convex optimization 
problem. A simple argument can show that the packet loss 
probability with the per-hop delay constraint is larger than 
that using the more flexible end-to-end constraint. 



B. Half-duplex Relay in Multihop Network 

Half-duplex relay is not a standard queue tandem model. 
However we can also derive a large deviation result for 
the tail probability for the end-to-end delay of a multihop 
network with half-duplex relays (proof in the Appendix): 

Theorem 5: For a stationary M/M/l queue tandem (with 
half-duplex relays), when the number of node N is large: 
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(17) 



From this theorem we conclude that the optimal ARQ alloca- 
tion problem with the end-to-end constraint and half-duplex 
relays can be formulated the same as that with full-duplex 
relays (14). 



IV. Numerical Examples 

Consider a MIMO relay network consists of a source, a 
relay, and a destination node. The relay is full-duplex. The 
number of antennas on each node is (Mi, M2, M3), Mi = 4, 
AI2 = 1, and M3 = 2, where the relay has a single antenna. 
Other parameters are: p = 20dB, k = 30, L = 8, and the 
multiplexing gain is r = 2. The base 10 logarithm of the 
cost function (7) is shown in Fig. 3. We have optimized the 
cost function with respect to L2 and k.2 so we can display it 
in three dimensions. Note that the surface is convex in the 
interior of the feasible region. The optimal L*, L\, k\ are 
shown in Fig. 4. Also note that as r increases to the maximum 
possible r = 4, the total number of ARQ rounds allocated 
L\ + L2 gradually increases to the upper bound L = 8 as k 
increases. 



Lindley's recursion has the solution: 

W* = max(4 in _i-4+i lT ,), i = 2,.--N-2, 



where the partial sum t\ v = A% and <r Up = £* =J S l k . 

Hence 
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From (20) we have r/ p 
and otherwise. Plug this into (23) we have 

Di = max^+i.! + Si - rj-\,„ - + D*" 1 ). (24) 

Hence the recursive relation if we move D l ~ Y to the left- 
hand-side: 



V. Conclusions and Future Work 

We have studied the diversity-multiplexing-delay tradeoff 
in multihop MIMO networks by considering an optimal 
ARQ allocation problem to minimize the probability of error, 
which consists of the ARQ error and the packet loss due to 
queueing delay. Our contribution is two-fold: we combine the 
DMDT analysis with queueing network theory, and we use 
the tail probability of random delay to find the probability 
of packet loss due to queueing delay. Numerical results 
show that optimal ARQ should equalize the performance 
of each link and avoid long service times that cause large 
queueing delay. Future work will investigate joint source- 
channel coding in multihop MIMO relay networks, extending 
the results of [5]. 

Appendix 

Proof of Theorem 5 

For node i, i = 1 • • • N, let the random variable S' l n denotes 
the service time required by the nth customer at the ith node 
(the number of ARQs used for the nth packet), and A? n be the 
inter arrival time of the nth packets (i.e., the time between the 
arrival of the nth and (n — l)th packages to this node). The 
waiting time of the nth packet at the ith node W£ satisfies 
Lindley's recursion (see [9]): 



Wl 



(W^_ 1 + S^_\-A l n )+, 2<z<7V_ 2 , (18) 



where (x) + = max(i, 0). The total time a message spent in 
a node is its waiting time plus its own service time, hence 



]j* =W l + S l 



(19) 



The arrival process to the (i + l)th node is the departure 
process from the ith node, which satisfies the recursion: 



At = Atr 1 + Dtr 1 - d 



71—1 ' 



2 < i < N -2. 



(20) 



with A l n a Poisson process with rate Also the waiting 
time at the source satisfies: 



Wl - {Wl_ x + + Sl_ x - Al)+. 



(21) 



A well-known result is that (see, e.g. [9]), if the arrival 
and service processes satisfy the stability condition, then the 
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Now from (23) we have D*, 1 = maxj^ ^(cr^. + 
Plug this in the above (25) we have 
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Do this inductively, we have 
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Ja <— < 3N - 1 =n 



If we also add D\ = W 7 \ + S^ to the above equation, after 
rearranging terms we have: 
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(26) 



Note that ct* ( ._ is independent of SJj . For long 

queue we can ignored the last four terms caused by edge 
effect (the source and end queue of the multihop relay 
network). By stationarity of the service process a\ ■ _ 1 + 
S) has the same distribution as a l n „■ . Then (26) 

3{i+l) Uji-J(i-l) v ' 

reduces to the case studied in [9] and we can borrow the 
large deviation argument therein to derive the exponent 9*. 
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Abstract — Tradeoff in diversity, multiplexing, and delay in 
multihop MIMO relay networks with ARQ is studied, where the 
random delay is caused by queueing and ARQ retransmission. 
This leads to an optimal ARQ allocation problem with per- 
hop delay or end-to-end delay constraint. The optimal ARQ 
allocation has to trade off between the ARQ error that the 
receiver fails to decode in the allocated maximum ARQ rounds 
and the packet loss due to queueing delay. These two probability 
of errors are characterized using the diversity-multiplexing- 
delay tradeoff (DMDT) (without queueing) and the tail proba- 
bility of random delay derived using large deviation techniques, 
respectively. Then the optimal ARQ allocation problem can be 
formulated as a convex optimization problem. We show that the 
optimal ARQ allocation should balance each link performance 
as well avoid significant queue delay, which is also demonstrated 
by numerical examples. 

I. Introduction 

In a multihop relaying system, each terminal receives 
the signal only from the previous terminal in the route 
and, hence, the relays are used for coverage extension. 
Multiple input-multiple output (MIMO) systems can provide 
increased data rates by creating multiple parallel channels and 
increasing diversity by robustness against channel variations. 
Another degree of freedom can be introduced by an automatic 
repeat request (ARQ) protocol for retransmissions. With the 
multihop ARQ protocol, the receiver at each hop feeds back 
to the transmitter a one-bit indicator on whether the message 
can be decoded or not. In case of a failure the transmitter 
sends additional parity bits until either successful reception 
or message expiration. The ARQ protocol provides improved 
reliability but also causes transmission delay of packets. Here 
we study a multihop MIMO relay system using the ARQ 
protocol. Our goal is to characterize the tradeoff in speed 
versus reliability for this system. 

The rate and reliability tradeoff for the point-to-point 
MIMO system, captured by the diversity-multiplexing trade- 
off (DMT), was introduced in [?]. Considering delay as the 
third dimension in this asymptotic analysis with infinite SNR, 
the diversity-multiplexing-delay tradeoff (DMDT) analysis 
for a point-to-point MIMO system with ARQ is studied in 
[?], and the DMDT curve is shown to be the scaled version 
of the corresponding DMT curve without ARQ. The DMDT 
in relay networks has received a lot of attention as well (see, 
e.g., [?].) In our recent work [?], we extended the point-to- 
point DMDT analysis to multihop MIMO systems with ARQ 
and proposed an ARQ protocol that achieves the optimal 
DMDT. 

The DMDT analysis assumes asymptotically infinite SNR. 
However, in the more realistic scenario of finite SNR, re- 
transmission is not a negligible event and hence the queueing 



delay has to be brought into the picture (see discussions in 
[?]). With finite SNR and queueing delay, the DMDT will be 
different from that under the infinite SNR assumption. The 
DMDT with queueing delay is studied in [?] and an optimal 
ARQ adapted to the instantaneous queue state for the point- 
to-point MIMO system is presented therein. 

In this work, we extend the study [?] of optimal ARQ 
assuming high but finite SNR and queueing delay in point- 
to-point MIMO systems to multihop MIMO networks. This 
work is also an extention our previous results in [?] to 
incorporate queueing delay. We use the same metric as that 
used in [?], which captures the probability of error caused by 
both ARQ error, and the packet loss due to queueing delay. 
The ARQ error is characterized by information outage prob- 
ability, which can be found through a diversity-multiplexing- 
delay tradeoff analysis [?], [?]. The packet loss is given 
by the limiting probability of the event that packet delay 
exceeds a deadline. Unlike the standard queuing models for 
networks (e.g., [?], [?]) where only the number of messages 
awaiting transmission is studied, here we also need to study 
the amount of time a message has to wait in the queue 
of each node. Our approach is slightly different from [?], 
where the optimal ARQ decision is adapted per packet; we 
study the queues after they enter the stable condition, and 
hence we use the stationary probability of a packet missing 
a deadline. An immediate tradeoff in the choice of ARQ 
round is: the larger the number of ARQ attempts we used 
for a link, the higher the diversity and multiplexing gain we 
can achieve, meaning a lower ARQ error. However, this is at 
a price of more packet missing deadline. Our goal is to find 
an optimal ARQ allocation that balances these two conflicting 
goals and equalizes performance of each hop to minimizes 
the probability of error. 

The remainder of this paper is organized as follows. 
Section II introduces system models and the ARQ protocol. 
Section III presents our formulation and main results. Nu- 
merical examples are shown in Section IV. Finally Section 
V concludes the paper. 

II. Models and Background 

A. Channel and ARQ Protocol Models 

Consider a multihop MIMO network consisting of 
nodes: with the source corresponding to i = 1, the destination 
corresponding to i = N, and i = 2, ■ ■ ■ , N — l corresponding 
to the intermediate relays, as shown in Fig. 1. Each node 
is equipped with Mj antennas. The packets enter the net- 
work from the source node, and exits from the destination 
node, forming an open queue. The network uses a multihop 



automatic repeat request (ARQ) protocol for retransmission. 
With the multihop ARQ protocol, in each hop, the receiver 
feeds back to the transmitter a one-bit indicator about whether 
the message can be decoded or not. In case of a failure 
the transmitter retransmits. Each channel block for the same 
message is called an ARQ round. We consider the fixed 
ARQ allocation, where each link i has a maximum of ARQ 
rounds Li, i = 1, • • • N — 1. The packet is discarded once 
the maximum round has been reached. The total number of 
ARQ rounds is limited to L: Xiila 1 — This fixed ARQ 
protocol has been studied in our recent paper [?]. 
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Fig. 1: Upper: relay network with direct link from source to 
destination. Lower: multihop MIMO relay network without direct 
link. 

Assume the packets are delay sensitive: the end-to-end 
transmission delay cannot exceed k. One strategy to achieve 
this goal is to set a deadline ki for each link i with 
Sil^ 1 k» — Once a packet delays more than ki it 
is removed from the queue. This per-hop delay constraint 
corresponds to the finite buffer at each node. Another strategy 
is to allow large per-hop delay while imposing an end-to-end 
delay constraint. 

Other assumptions we have made for the channel models 
are 

(i) The channel between the ith and (i + l)th nodes is given 
by: 



Y h i - J -^H u Xu + W u , l<l< Li. (1) 

The message is encoded by a space-time encoder into 
a sequence of L matrices £ C MiXT ,: I = 

1, • ■ • , L}, where T is the block length, and Yij G 
C M,+lxT , i = 1, ■ • • , N — 1, is the received signal at 
the (i + l)th node, in the Zth ARQ round. The rate 
of the space-time code is R. Channels are assumed 
to be frequency non-selective, block Rayleigh fading 
and independent of each other, i.e., the entries of the 
channel matrices Hij G C M '+ lxMi are independent 
and identically distributed (i.i.d.) complex Gaussian with 
zero mean and unit variance. The additive noise terms 
Wij are also i.i.d. complex Gaussian with zero mean 
and unit variance. The forward links and ARQ feedback 
links only exist between neighboring nodes. 

(ii) We consider both the full-duplex and half-duplex relays 
(see, e.g., [?]) where the relays can or cannot transmit 
and receive at the same time, respectively, as shown 
in Fig. 2. Assume the relays use a decode-and-forward 
protocol (see, e.g., [?]). 

(iii) We assume a short-term power constraint at each node 
for each block code. Hence we do not consider power 
control. 

(iv) We consider both the long-term static channel, where 
Hi.i = Hi for all I, i.e. the channel state remains 



constant during all the ARQ rounds, and independent for 
different i. Our results can be extended to the the short- 
term static channel using the DMDT analysis given in 
[?]. 
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Fig. 2: Left: full duplex multihop relay network. Right: half duplex 
relay multihop MIMO relay network. 




Fig. 3: The logarithm of the cost function (7) for the (4, 1, 2) 
multihop MIMO relay networks. SNR is 20 dB. 
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Fig. 4: Allocations of optimal ARQ: L\, L* 2 , LI+L2, k{, in a (4, 1, 
2) MIMO relay network. SNR is 20 dB. (The optimal fcj = fc-fci-) 

B. Queueing Network Model 

We use an M/M/l queue tandem to model the multihop 
relay networks. The packets arrive at the source as a Poisson 
process with mean interarrival time fi, (i.e., the time between 
the arrival of the nth packet and (n — l)th packet.) We 
model the random service times for the nth packet at the zth 
node as independent and identically exponentially distributed 
with mean Li that is also independent of the inter-arrival 
time. This approximation leads to a tractable solution that 
characterizes the qualitative behavior of the multihop network 
queue. Node m has a finite buffer size. The packets enter into 
the buffer and are first-come-first-served (FCFS). Assume 
/i > Li so that the queues are stable, i.e., the waiting time 
at a node does not go to infinity as time goes on. 



Burke's theorem (see, e.g., [?]) says that the packets depart 
from the source and arrive at each relay as a Poisson process 
with rate p m /p, where p$ is the probability that a packet 
can reach the mth node. With high SNR, the packet reaches 
the subsequent relays with high probability: pj « 1 (the 
probability of a packet dropping is small because it uses up 
the maximum ARQ round.) Hence all nodes have packets 
arrive as a Poisson process with mean inter-arrival time p. 

C. Throughput 

Denote by 6 the size of the information messages in bits, 
B[t] the number of bits removed from transmission buffer 
at the source at time slot t. Define a renewal event as the 
event that the transmitted message leaves the source and 
eventually is received by the destination node possibly after 
one or more ARQ retransmissions. We assume that under 
full-duplex relays the transmitter cannot send a new message 
until the previous message has been decoded by the relay at 
which point the relay can begin transmission over the next 
hop (Fig. 2a.) Under half duplex relays we assume transmitter 
cannot send a new message until the relay to the next hop 
completes its transmission (Fig. 2b.) 

The number of bits B transmitted in each renewal event, 
for full-duplexing B = (N — 1)6, and for half-duplexing 
B = (N— 1)6/2 when N is odd, and B = Nb/2 when N is 
even. The long-term average throughput of the ARQ protocol 
is defined as the transmitted bits per channel use (PCU) [?], 
which can be found using renewal theory [?]: 
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R, Full duplex; 

f , Half duplex, N is odd; (2) 

R (i + jjy) , Half duplex, TV is even. 

where r is the average duration from the time a packet arrives 
at the source until it reaches the destination node, and = 
denotes asymptotic equality. A similar argument as in [?] 
shows that E(t) = (N - 1)T for high SNR. 

D. Diversity-Multiplexing-Delay Tradeoff 

The probability of error P e in the transmission has two 
sources: from the ARQ error: the packet is dropped because 
the receiver fails to decode the message within the allocated 
number of ARQ rounds, denoted as P ARQ , and the probability 
that a message misses its deadline at any node due to large 
queueing delay, denoted as P Que „ t . We will give P e for various 
ARQ relay networks. 

Following the framework of [?], we assume the size of 
information messages b(p) depends on the operating signal- 
to-noise ratio (SNR) p, and a family of space time codes 
{C p } with block rate R(p) = b{p)/T = rlogp. We use the 
effective ARQ multiplexing gain and the ARQ diversity gain 
[?] 



lim 



d = — lim 

log p p^oo log p 

We cannot assume infinite SNR because otherwise the queue- 
ing delay will be zero, as pointed out in [?]. However we 
assume high SNR to use the DMDT results in our subsequent 
analysis. 
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III. Diversity, multiplexing, and delay tradeoff 

VIA OPTIMAL ARQ ROUND ALLOCATION 
A. Full-Duplex Relay in Multihop Relay Network 

1) Per-Hop Delay Constraint: The probability of error 
depends on the ARQ window length allocation Li, deadline 
constraint h- L , multiplexing rate r, and SNR p. For a given r 
and p, we have 



P e ({Li},{ki}\ P ,r) 



N-l 



P ARQ ( P , {U}) + P^(Di > h). 



(4) 



Here Di denotes the random delay at the ith link when the 
queue is stationary. This P e expression is similar to that 
given by Equation (33) of [?]. Our goal is to allocate per-hop 
ARQ round {Li} and delay constraint {ki} to minimize the 
probability of error P e . 

For the long-term static channel, using the DMDT analysis 
results [?1 we have: 
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P ARQ (p, {U}) = J2 P f '^' 



(5) 



Here fi(r) is the diversity-multiplexing tradeoff (DMT) for 
a point-to-point MIMO system formed by nodes i and i + 
Assuming sufficient long block lengths, fi(r) is given by 
Theorem 2 in [?] quoted in the following: 

Theorem 1: [?] For sufficiently long block lengths, the 
diversity-multiplexing tradeoff (DMT) f(r) for a MIMO 
system with M t transmit and M r receive antennas is given 
by the piece-wise linear function connecting the points 
(r, (M t - r)(M r - r)), for r = 0, • • • , mm(M t , M r ). 

Denote the amount of time spent in the ith node by 
the ?ith message as D z n . The probability of packet loss 
Poueuc(-Di > h) can be found as the limiting distribution of 
lim^oo P(D i n > h) (adapted from Theorem 7.4.1 of [?]): 

Lemma 2: The limiting distribution of the event that the 
delay at node i exceeds its deadline ki, for M/M/l queue 
models, is given by: 



P<t~*(Di > ki) 
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(6) 



Here the difference in the service rate and packet arrival rate 
j- — — > and utility factor both indicate how "busy" 
node m is. 

Using the above results, (4) can be written as 



P e {{L i },{k i }\p,r) = 
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(7) 



Note that the queueing delay message loss error probability is 
decreasing in Lj, and the ARQ error probability is increasing 
in L^ Hence an optimal ARQ rounds allocation at each node 
Li should trade off these two terms. Also, the optimal ARQ 
allocation should also equalize the performance of each link, 
as the weakest link determines the system performance [?]. 

Hence the optimal ARQ allocation can be formulated as 
the following optimization problem: 



„ mm Pe({Li},{h}\p,r) 
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where 

f T^L^L, 
A = < l<L,<p, i = !,••• ,N — I 

I EtT 1 fa < *• 

The following lemma (proof omitted due to the space limit) 
shows that the total transmission distortion function (15) is 
convex in the interior of A. 

Lemma 3: The transmission distortion function (15) is 
convex jointly in Li and k{ in the convex set 
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Lemma 3 says that except for the "corners" of A the cost 
function is convex. However these "corners" have higher 
probability of error: fej and Li take extreme values and hence 
one link may have a longer queueing delay then the others. 
So we only need to search the interior of A where the cost 
function is convex. 

To gain some insights into where the optimal solution 
resides in the feasible domain for the above problem, we 
present a marginal cost interpretation. Note that the proba- 
bility of error can be decomposed as a sum of probability 
of error on the zth link. The optimal ARQ rounds allocated 
on this link should equalize the "marginal cost" of the ARQ 
error and the packet loss due to queueing delay. For node i, 
with fixed fcj, the marginal costs (partial differentials) of the 
ARQ error probability, and the packet loss probability due to 
queueing delay, with respect to Li are given by 
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dLi p. 

Note that f! < 0. The optimal solution equalizes these 
two marginal costs by choosing Lj G Note that 

these marginal cost functions are monotone in Li, hence the 
equalizing L* exists and 1 < L\ < p if the following two 
conditions are true for L\ = 1 and Li = p: 



(i): 
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These conditions involve nonlinear inequalities involving 
p, p, r, Mi and Mi+i, which defines the case when the 
optimal solution is in the interior of A. Analyzing these 
conditions reveals that these conditions tend to satisfy at 
lower multiplexing gain r, small Mj or Mj+i, small foj, and 
larger p (light traffic). Note that with high SNR condition 
(ii) is always true for moderate k values. When (i) and (ii) 
are violated, which means one error dominates the other, 
then the optimal solution lies at the boundary of A. With 
the total ARQ rounds constraint in (8), using the Lagrangian 
multiplier an argument similar to above still holds. 



2) End-to-End Delay constraint: When the buffer per 
node is large enough a per hop delay constraint is not needed, 
and we can instead impose an end-to-end delay constraint. 
The exact expression for the tail probability of the end-to- 
end delay is intractable. However a large deviation result is 
available. The following theorem can be derived using the 
main theorem in [?]: 

Theorem 4: For a stationary M/M/l queue tandem (with 
full-duplex relays): 
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where 



This theorem says that the bottleneck of the queueing net- 
work is the link with longest mean service time Hence 
the optimal ARQ round allocation problem can be formulated 

as: 



min P e ({Li},{ki}\p,r) 
{Li}eB 



(14) 



where 



P e ({£,}, {fci}|p,r) 
= P ARQ (p, {Li}) + P Q11CUC ( > k ) , 
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For high SNR, this can be shown to be a convex optimization 
problem. 

A simple argument can show that the packet loss probabil- 
ity with the per-hop delay constraint is larger than that using 
the more flexible end-to-end constraint. 

B. Half-duplex Relay in Multihop Network 

We can also find the tail probability for the end-to-end 
delay of a multihop network with half-duplex relays using 
large deviation techniques (see proof in the Appendix). Thus 
yields 

Theorem 5: For a stationary M/M/l queue tandem (with 
half-duplex relays), when the number of node N is large: 
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From this theorem we conclude that the optimal ARQ alloca- 
tion problem with the end-to-end constraint and half-duplex 
relays can be formulated the same as that with full-duplex 
relays (8) and (14). 

IV. Numerical Examples 

Consider a MIMO relay network consists of a source, a 
relay, and a destination node. The relay is full-duplex. The 
number of antennas on each node is (Mi, M2, M3), Mi = 4, 
M2 = 1, and M3 = 2, where the relay has a single antenna. 
Other parameters are: p = 20dB, k — 30, L = 8, and the 



multiplexing gain is r = 2. The base 10 logarithm of the 
cost function (7) is shown in Fig. 3. We have optimized 
the distortion function with respect to L2 and k% so we can 
display it in three dimensions. Note that the surface is convex 
in the interior of the feasible region. The optimal L\, L\, k\ 
are shown in Fig. 4. Also note that as r increases to the 
maximum possible r = 4, the total number of ARQ rounds 
used gradually increases to the upper bound L = 8 as k 
increases. 

V. Conclusions and Future Work 

We have studied the diversity-multiplexing-delay tradeoff 
in multihop MIMO networks by considering an optimal 
ARQ allocation problem to minimize the probability of error, 
which consists of the ARQ error and the packet loss due to 
queueing delay. Our contribution is two-fold: we combine the 
DMDT analysis with queueing network theory, and we use 
the tail probability of random delay to find the probability 
of packet loss due to queueing delay. Numerical results 
show that optimal ARQ should equalize the performance 
of each link and avoid long service times that cause large 
queueing delay. Future work will investigate joint source- 
channel coding in multihop MIMO relay networks, extending 
the results of [?]. 

Appendix 

Proof of Theorem 5 

For node i, i = 1 ■ • • N, let the random variable 5* denotes 
the service time required by the nth customer at the ith node 
(the number of ARQs used for the nth packet), and A % n be the 
inter arrival time of the nth packets (i.e., the time between the 
arrival of the nth and (n — l)th packages to this node). The 
waiting time of the nth packet at the ith node W„ satisfies 
Lindley's recursion (see [?]): 



From (20) we have r. 



l, P 



•l, P 



D 



i-l 



D\z\ for / <p + l, 



and otherwise. Plug this into (23) we have 
Di = max(<7*+;_i + S l n - r^ l n - D^ 1 + D) 1 1. (24) 

Hence the recursive relation if we move D^ 1 to the left- 
hand-side: 

Di + D*- 1 = ma X (a}+^ 1 + - r£ 1>n + D) 1 i. (25) 
Now from (23) we have D l -~ l = 
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Plug this in the above (25) we have 
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Do this inductively, we have 
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If we also add D\ = + to the above equation, after 
rearranging terms we have: 
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(26) 



Note that a\ 4 , is independent of S] . For long 
queue we can ignored the last four terms caused by edge 
effect (the source and end queue of the multihop relay 
network). By stationarity of the service process a 
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SI has the same distribution as a % n ■ ■ . Then (26) 
reduces to the case studied in [?] and we can borrow the 
large deviation argument therein to derive the exponent 6*. 
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where (x) + = m&x(x, 0). The total time a message spent in 
a node is its waiting time plus its own service time, hence 
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The arrival process to the (i + l)th node is the departure 
process from the ith node, which satisfies the recursion: 



A' 1 — A 1 
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(20) 



with A l n a Poisson process with rate 1/fj,. Also the waiting 
time at the source satisfies: 
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A well-known result is that (see, e.g. [?]), if the arrival 
and service processes satisfy the stability condition, then the 
Lindley's recursion has the solution: 
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where the partial sum t\ p = YX=i A k and a i,P = Yl=i S l k - 
Hence 



Dm — m 5 x ( (J i«,n-i + T ji+1. 



i = 2,---M-2. (23) 



